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Recently, dissipationless spin current has been predicted for the p-doped semiconductors with spin- 
orbit coupling. Here we investigate the effect of spherical symmetry breaking on the dissipationless 
spin current, and obtain values of the intrinsic spin Hall conducitivity for realistic semiconductor 
band structures with cubic symmetry. 
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Spintronics is a new field of science and technology 
which aims to manipulate the spin of the electron for 
building functional logic and storage devices Q. The 
creation, manipulation and transport of spin currents is 
a central challenge in this field. Recently, Murakami, 
Nagaosa and Zhang|2| found a basic law of spintronics, 
which relates the spin current and the electric field by 
the response equation 

j] = (TseijkEk (1) 

where j] is the current of the i-th component of the spin 
along the direction j and Cijk is the totally antisymmet- 
ric tensor in three dimensions. This effect arises because 
of the spin-orbit coupling in the valence band of con- 
ventional semiconductors such as GaAs and Ge. Sinova 
et al0 also found a similar effect in the electron doped 
conduction band. The transport equation Q is similar 
to Ohm's law in electronics. However, unlike the Ohm's 
law, this new law describes a purely dissipationless spin 
current, in the sense that Eq. ^ is invariant under the 
time reversal and the intrinsic part of as does not depend 
on impurity scattering. These effects have been further 
discussed in the recent literature 0, 0, 0, 13. 0. • 

Fundamental to the proposal of Murakami, Nagaosa 
and Zhang Q is the spin-orbit coupling that exists in the 
Luttinger effective-mass model in the degenerate valence 
bands: 

H^^(i7i + ll2)k'-2j,{k-Sf^ (2) 

In this "isotropic," or spherically symmetric model, the 
helicity A — k ■ S is a good quantum number of the 
isotropic Luttinger hamiltonian above, and it labels the 
two doubly degenerate Kramers bands corresponding to 
the heavy holes A = ±| and light holes A = ±i . The 
spin current effect can be intuitively understood as a con- 
sequence of the conservation of total angular momentum: 
J = ?ix X k-|- S. The spin current fiows in such a way that 
the change of the orbital angular momentum L = ftx x k 
exactly cancels the change of the spin angular momentum 



S. When an electric field is applied on the arbitrary z 
axis, the z component of J is conserved. The topological 
nature of the spin current is manifest in th e gauge-field 
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Si 


4.22 


0.39 


1.44 


0.248 


Ge 


13.35 


4.25 


5.69 


0.108 


GaAs 


4.15 


1.01 


1.75 


0.114 


InSb 


35.08 


15.64 


16.91 


0.036 


InAs 


19.67 


8.37 


9.29 


0.047 


GaP 


4.20 


0.98 


1.66 


0.162 



TABLE 1: Valence-band parameters for some common ma- 
terials [Tol . Following [Tl| we define S = (73 — 72)/7i as a 
measure of the anisotropy. 

formulation of , where the spin conductance is defined 
in terms of a linear combination of the components of a 
gauge field, Gij — A(A^ — 13 /4:)eijiki/k^, clearly reflect- 
ing a monopole structure in k space. The singularity 
at k — > exemplifies the confluence of the Kramers 
doublets at the F point where the band becomes 4 fold 
degenerate, but the flux of the gauge field through a 2 
dimensional surface in k space is constant and set by the 
helicity eigenvalue. 

The picture presented above is valid as long as the 
Hamiltonian is isotropic, that is to say it has spherical 
symmetry. In the real materials in which the dissipation- 
less spin-current is predicted '^l , all of which are charac- 
terized by a large anisotropy (see Table QJ, the angular 
momentum J and the helicity X — k-S are no longer good 
quantum numbers. It is therefore vital to ask whether the 
topological spin current is preserved in materials which 
are not rotationally invariant. In this paper, we investi- 
gate the effect of the spherical symmetry breaking on the 
dissipationless spin current, and calculate the values of 
the intrinsic spin Hall conductivity for anisotropic band 
structure parameters. 

The most general Hamiltonian which respects time- 
reversal and cubic symmetries was derived by Lut- 
tinger [1^ : 
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Ho = T^(7i + ^72)fc' - —{klSl + klSl + klSl) - 2^[{fc,, fcj{5,, Sy} + {fc^, fcj{^j„ 5J + {K,K}{S,, S,}] 
2m 2 m " " m 

(3) 



where we define {A, B} = ^{AB + BA), and ^ kl + 
ky + k'^. The parameters, 71,72, and 73, are material- 
dependent. In the special case of 72 — 73 (which we call 
isotropic), the last two terms simply combine to yield 

In real materials, however, the values of 72 and 73 are 
very different. Table 1 lists the values of these parameters 
in some important materials. The anisotropy, character- 
ized by the parameter S = (73 — 72)/7i, is relevant and 
substantial for all the materials, and especially relevant 
for Si. In order to understand the dissipationless spin- 
current generated in these real materials, including its 
dependence on the orientation of the field and current 
with respect to the crystal axes, we must consider the 
full anisotropic Hamiltonian, Eq. 

When 72 7^ 73, the Hamiltonian is no longer isotropic 
and the helicity is not a good quantum number. However, 
the energy spectrum of the Hamiltonian retains the same 
structure as in the isotropic case, albeit with a different 
dispersion relation. After diagonalizing the Hamiltonian, 
we obtain two doubly degenerate energy levels, which we 
call light and heavy holes in analogy with the isotropic 
case: 



E{k) ^^-fik^±—d{k) 
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(4) 



Following Ref. jSj, we can expand the spin-dependent 
terms in the anisotropic Luttinger Hamiltonian in terms 
of a Clifford algebra of Dirac F matrices {F°,F''} = 
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Ho = e(k) + ^dj" 



(5) 



di — —VSkzky, d2 — ' 

2 73 



yf^kx ky ^ 



d, = -\':il{2kl-kl 
273 
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with dada = . Whereas in the isotropic Luttinger 
model the matrix used to diagonalize the Hamiltonian 
belongs to the 5*0(3) group of rotations in k space 2], in 
the anisotropic materials the matrix that diagonalizes the 



anisotropic Hamiltonian belongs to the 50(5) rotations 
in da space. The 50(5) Clifford algebra representation of 
the Hamiltonian |SJl naturally unifies both the isotropic 
and the anisotropic Luttinger model on the same footing. 
Since this form of the Hamiltonian depends on k only 
through the five dimensional vector da, a large part of 
the results in is directly applicable to the anisotropic 
case. In this sense, the S'0(5) Clifford algebra formal- 
ism shows its full power in the anisotropic case studied 
here. The projection operators onto the two-dimensional 
subspace of states of the heavy-hole (HH) and light-hole 
(LH) bands read: 

P^ = ^(l + dar") , P^ = i(l-d;F'^) (7) 

For finite k, the Hamiltonian maintains the 50(4) sym- 
metry observed in 5] . This symmetry reflects the degen- 
eracy of the two Kramers doublets at each value of k, 
corresponding to the doubly-degenerate HH and the LH 
bands. Each of the bands has an SU (2) symmetry, which 
we may denote by SU{2)hh and SU{2)lh- Therefore, 
the total symmetry is SU{2)hh ^ SU{2)lh = 50(4). At 
the r point, k = 0, there is an enhanced 50(5) symme- 
try. 

The symmetry generators read: 
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ab 



dhdrV" 



pL^dh pL 



-^cb 



pHYabpH 



(8) 



where F"'' = -|[r",F*'] and i?o] = trivially since 
the Hamiltonian is diagonal in the HH and LH bands. 
The spin operators 5* are related to the F"'' matrices 
through the tensor ry*^, whose entries were given in j^: 
5* = Vab^'^''- The concept of a conserved spin cur- 
rent is still valid in anisotropic materials, since the pro- 
jected spin is a constant of motion in virtue of its be- 
ing a linear combination of the symmetry generators, 
5|^^ = ril^p"'' = P^S^P^ + P^S'P". We can therefore 

define the conserved spin-current as = ^ | ^ , 5|^j | . 
Note that the richer anisotropic Luttinger Hamiltonian 
yields a very similar structure to the isotropic one when 
cast in 50(4) language. 

Although the concept of helicity A = /ci5* is not valid 
in anisotropic materials, we can define a corresponding 
conserved helicity, Xnew, as: 



At, 



ki^lc) — A 



2hril,dhd,T-'' = 



pH^pH 



(9) 
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Since it is a linear combination of the symmetry gener- 
ators of Hq {Xnew = fci'S'(c) = kiilabP"''')^ clcar that 
[iJ, Xnew] — 0. In the isotropic hmit, A„e„, = A, as can 
be seen using the identities [A, P^] = [A, P^] = 0, vahd 
in the isotropic case. 

The recent work of Ref.j^ shows that the Kubo for- 
mula for the conserved spin current response can be ex- 
pressed purely in terms of a geometric quantity 



/-^ f^abj^ab 

"^ij — '^ij ^ — ^abcde 



ddd ddf. 

73^ea6c.e4^^r (10) 



which describes the mapping from the 3D k vector space 
to the 5D d vector space. This results also includes a 
quantum correction to the semiclassical result of Ref . . 
We shall apply this formula to the anisotropic case here. 
However, there is one essential difference. Whereas in 
the isotropic case, the field strength can be brought, 
through a proper choice of gauge, to the diagonal form 
Gij = A(A^ — 13 /A)eijiki/k^, in the anisotropic case 
this is impossible. Non-abelian field strength are, in 
general, gauge- variant. However, there is a fundamen- 
tal difference between fields that can be diagonalized 
through a gauge transformation and fields for which this 
is not possible. The former are ultimately abelian in 
nature, whereas the latter are truly non-abelian. The 
non-diagonal gauge field which describes evolution in 
anisotropic materials reflects the richer structure of the 
anisotropic Luttinger hamiltonian. 

We can express the field strength in terms of the (un- 
projected) spin degrees of freedom if we first note that 
the ten SO{5) generators F"'' decompose into the 3 spin 
matrices and the seven cubic, symmetric and traceless 
combinations of the spin operators of the form S^S^ S'^, 
namely: 



A^ = {Sx,{Syr-iSzn 

A^^{Sy,{S,y'-{Sx)^} 
A^^iSzASxr-iSy)'} 

A. — SxSySz 



(11) 



SzSySx 



Then we can write: 



G^, = ^e.uMV^A'^ + U,S' 



I = 1, 



,3,^ = 1, 



-J 
(12) 



Where 



Ui - i^[(13 



28^)fcf 
~3^A:a,(/Cy 
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V, = ^3fkz{kl 
V7 — Xlkxkykz 



-28^)fc2fc, 

A^)Pki 
73 ' ' 



73 ' 



(13) 



When cast in the 50(4) language, the expression for 
the spin conductance in anisotropic materials has the 
same form as in the spherical model: 



Vh 



Y,{nLiy)-nHm\r,\,G'i^ (14) 



where ul = npi^L) and nu = np{€H) are the Fermi 
functions of the LH and HH bands. This expression can 
be put into the following elegant form: 
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where we see that the first term in brackets vanishes in 
the isotropic case. The I index specifies the direction 
of the spin orientation, and it is not summed over on 
the right hand side of Eq. (|15() . It is now obvious that 
the only components of cr'^- surviving after summing the 
contributions from the whole Fermi surface, are those 
for which i ^ j ^ I. Indeed, upon integration over k, 
ct'j becomes proportional to Ci,^ just as it should for in 
crystals with cubic symmetry |l3l | . 

Our result for the spin-current can thus be put in the 
form: 



CTs = -^"-^^^'S'(7i,72,73), 



(16) 



Z = 1,..,3 



where the material-specific coefficient, S, is independent 
of the Fermi energy, and is of the order ~ 0.05 for 
most materials (see Table Hi)) . The Cs ''^ n^^^ scaling 
is the hallmark of the dissipationless spin current, and 
has been proposed as a means to distinguish from other 
extrinsic effects 0, 0- To compare the spin-conductance 
in different materials, we separate the dependence on 
the total carrier density, which for the anisotropic Lut- 
tinger model depends on the band-parameters: n — 
(Omp f \ 1 1 i 1 t^^fc 

y^iiLt-F) 3 J [(.^i_273d/fc2)3/2 -r (^j+2-),3d/fc2)3/2j IS^P"- 

Using this relation, we can find ep as a function of n, 
and use it to define the anisotropic Fermi distribution 
functions, nL,_f/(k) = Q{kp'^ — k). We have calculated 
CTs for band parameters corresponding to a selection of 
real materials, as well as for band parameters correspond- 
ing to isotropic materials with the same values of 71 and 
11 = The results, listed in Table HH show that 

the non-zero anisotropy leads to a decrease in the spin- 
conductivity of as much as 30% (for Si), although the re- 
duction in materials with smaller anisotropy is typically 
only - 5%. 

To illustrate the systematic dependence of spin- 
conductance on anisotropy, we plot cTg as a function of 
5 = (73 — 72)771 with 71 and ji = (673 -I- 472)/5 held 
fixed at the values corresponding to Si, GaAs and InSb 
(Figure^]) . The spin-conductance at fixed carrier concen- 
tration is maximum at 5 = 0, whereas all real materials 
have 5 > 0. This observation should guide the selection of 
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5(71,72,73) 


as cm"^) 




Si 


0.028 


14.60 


21.10 


Ge 


0.063 


32.79 


34.31 


GaAs 


0.062 


32.64 


34.33 


InSb 


0.083 


43.63 


44.67 


InAs 


0.079 


41.61 


42.55 


GaP 


0.051 


26.50 


29.17 



TABLE II: We list the material-dependent coefficients of 
the spin-conductivity for values of 71,72,73 corresponding 
to common semiconductors. Also given are the actual spin- 
conductivities at n = lO^^cm""^ for both the real anisotropic 
materials, and their spherical approximations {5 = 0). 
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FIG. 1: Spin-conductivity is plotted as a function of the 
anisotropy, parameterized hy 5 = (72 — 73)/7i, with fj. = 
(673 -I- 472)75 and n = lO^^cm""^ held fixed at values corre- 
sponding to Si (bottom curve), GaAs and InSb (top curve). 
The circles indicate the real values of the parameters in Si, 
GaAs and InSb. 

materials with relatively low anisotropy for spin-injection 
devices and other applications where a stron spin-current 
is desired. Finally, the variation of as with carrier con- 
centration, 71, is shown in Figure |2| 
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FIG. 2: The dependence of spin-conductivity on carrier den- 
sity, n, is plotted using the band parameters of Si (bottom 
curve), GaAs and InSb (top curve). 
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